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MORPHOLOGICAL STUDIES OF RISING EQUATORIAL SPREAD F BUBBLES 


■ 


I. INTRODUCTION 

Recently, experimental ( Kelley et al. , 1976; Woodman and La Hoz , 
1976; and McClure et al. , 1977) and numerical simulation studies 
( Scannapieco and Ossakow , 1976) of the nighttime equatorial F region have 
given evidence for rising bubbles (plasma density depletions). The ver¬ 
tical speed of these bubbles has been observed ( McClure et al. , 1977) to 
range from very low (several meters/sec and even zero) to very high 
(MOO m/sec or more) velocities. Hudson [1977] has attempted to 
describe these rising bubbles within the context of the collisionless 
Rayleigh-Taylor mechanism using techniques in two dimensions that 
Chaturvedi and Kaw (1975 a,b) used in one dimension. Ott [1977] has 
described the rise of equatorial Spread F bubbles, in both the col 1i- 
sional and collisionless Rayleigh-Taylor regimes, by noting the analogy 
with bubbles in fluids. Szuszczewicz [1977] has discussed the chemistry 

and transport of equatorial Spread F ionospheric holes. 

It has been pointed out ( Scannapieco and Ossakow , 1976; and 
Scannapieco et al. , 1976) that there is a considerable analogy between 
equatorial Spread F (in the collision-dominated regime) and barium cloud 
studies. Indeed barium clouds can be viewed as Pedersen conductivity or 
plasma density enhancements while natural equatorial Spread F bubbles 
(or artificially created ionospheric holes, e.g., see Mendillo et al. , 

Note: Manuscript submitted November 7, 1977. 
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1975; and Bernhardt and da Rosa, 1977) can be viewed as Pedersen con¬ 


ductivity or plasma density depletions. In the present paper, we wish 

to discuss equatorial Spread F bubble shapes and vertical rise rates 

within the context of the collisional Rayleigh-Taylor regime (as invoked 
by Scannapieco and Ossakow , 1976), making use of the electrodynamic 

analogy (i.e., the electrical properties of the bubbles) with barium 

clouds. In section II we present the general theory, with subsections 

dwelling on linear and nonlinear bubble rise velocity models. Section 

III deals with the discussion and summary. 

II. THEORY 


The equations governing equatorial Spread F (within the context of 
the collisional Rayleigh-Taylor instability) are (see Scannapieco and 
Ossakow, 1976) 


I +v . 


T<?ax£ 


V tjj x b) 
1 


= -v R (n-n ) 


( 1 ) 


V i ( v in nV i^) = e 9- ' V v in n ^ + x ^ ’ v . n 

where c is the speed of light, B the magnetic field (taken to be con¬ 
stant), e the electronic charge, M the mass of the ions, £ acceleration 

A 

of gravity, b = B/|B| (pointing North), the subscript i means perpendi- 

A / 

cular to b (Note: the geometry we will be employing for equatorial 
Spread F is such that x is upward, y is horizontal and eastward and z 
is northward and along the earth's magnetic field), is the recombina¬ 
tion rate, n the ion number density (or equivalently a field line inte¬ 
grated number density in a slab geometry), E = -V^, where ^ is the 
electrostatic potential (in this model there is no ambient, i.e., zero 
order, horizontal electric field), v. the ion-neutral collision 
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frequency (a function of altitude) and n the equilibrium ionospheric 
electron (ion) number density. 

If we neglect recombination effects, neglect the variation of v. 
with altitude, and set = <f> + 4> (where from (1) and (2) V<p Q = ^), 
then we have 


fjr - | (V a 4 > x b) • V A n = 0 ( 3 ) 

V • (nV <f>) = - g x b • V n (4) 

i i v. e ■ 2 - i ' ' 

^ n 

(Note: The same potential equation (4) can be obtained by keeping terms 
of o(l) compared with v. /fi, where SI = eB/Mc, in eq. (2)). The simplest 
set of equations describing F region barium cloud morphology and stria- 
tions (for a one level model and in the E x B drift frame) is (see 
Perkins et al. , 1973 and Zabusky et al. , 1973) 


9 N c 

at “ B 



A 

x b) • V N = 0 

i 



V (N V <{>) = E • V N 

1 ' i T ' —0 i 



where is taken to be the ambient electric field perpendicular to B 

(E = E - V,<f>) and N is the magnetic field line integrated Pedersen 
—0 1 
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conductivity or equivalently the magnetic field line integrated plasma 
density. (Note: In these equations v. is taken to be a constant and 
so represents a rotated geometry, i.e., the barium cloud equations are 
more applicable to high latitudes where B is vertical, V N and E are 

— i -o 

horizontal.) If gravitational effects were also included in (6) then 

A 

E* e E + (ft/v,- )(M/e) q x b would replace E. on the right hand side of 

U I II U 

(6) (see Perkins et al. , 1973). The two sets of equations (3) and (4), 

A 

and (5) and (6) are essentially the same with (ft/v. )(M/e) g x b playing 
the role of an electric field in the potential eg. (4) (recall also the 
definition of E_*). 

Let us now discuss some depletion and enhancement cloud morphology. 
The separate coupled sets of equations (3) and (4) and (5) and (6) have 
no steady state solutions if the cloud, as described by n or N, has a 
finite size (with the gradient of n or N finite) in two dimensions per¬ 
pendicular (x,y) to the magnetic field (see Perkins et al. , 1973 and 
Dungey , 1958). Another result concerns the initial deformation of a 
cylindrical plasma density or Pedersen conductivity enhancement (barium 
cloud) and a cylindrical plasma density or Pedersen conductivity deple¬ 
tion (equatorial Spread F bubbles or artificial lycreated ionospheric 
holes). The situation is illustrated in Figure 1. Here in this context 
eq. (4) or (6) is simply the problem of cylindrical dielectric immersed 
in a uniform electric field (where in the equatorial Spread F situation 

A. 

(ft/v. )(M/e) £ x b plays the role of an equivalent electric field . In 
the case (A) of the density enhancement (e.g., barium cloud) the induced 
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electric field is opposed to the uniform electric field, so that E/, 

the total electric field inside the enhancement is less than E . In the 

-o 

case (B) of the density depletion (e.g., equatorial Spread F bubble or 
artificially created "hole") the induced electric field is in the direc¬ 
tion of the uniform electric field, E . so that E', the total electric 

—o — 

field inside the depletion is greater than E^. (Note: the formulas for 
E_' are approximate and specific models will be discussed later on in the 
text.) It should be pointed out that in Fig. la the circles representing 
the isodensity contours are such that the enhancement and depletion are 
maximum in the center and fall off as one moves away from the center. 

In this way the induced electric field goes to zero outside of the 
enhancement or depletion (Note: for the enhancement the induced electric 
field is — 6NE /(N + 6N), where N q is the ambient ionospheric density 
and 6N is the enhanced density, while for the depletion the induced 
electric field is 6nE /(n - 6n) where fin is the magnitude of the deple¬ 

tion and n is the ambient density). In case (A) the x drift of 
the center of the enhancement counteracts the general E' x B drift and 
the center of the enhancement (cloud) has the slowest drift. The back¬ 
side (side opposite to E x B motion) of the enhancement catches up with 

0 

the center and a steep density gradient forms there (depicted in Fig. lb). 
Similarly, the frontside elongates. The area enclosed by any isodensity 
contour remains constant during the motion and the maximum value of |N| 
is a constant of the motion (incompressibility). 

In case (B), the E 1 x B drift of the center of the depletion (bubble) 
acts in concert with the general E' x B drift and the center of the 
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depletion has the fastest drift. The center now catches up to the front¬ 
side and a steep density gradient forms there (depicted in Fig. lb). In 
this case the so-called backside (using barium cloud terminology) elongates. 
In the context of eq. (3) and (4) in Fig. 1, the area enclosed by any 
isodensity contour remains constant during the motion and the maximum 
value of |n| is a constant (incompressibility). In both the enhancement 
and depletion cases the density piles up in the direction formed by the 
induced electric field crossed with the magnetic field. In both cases 
where these steep gradients form instabilities should be generated. In 
some sense the density depletion can be likened to an "antibarium cloud." 
Within the context of Fig. 1 we have spoken about depletions in general 
and the general results apply to natural Spread F bubbles or artificially 
created ionospheric holes (using H-O, Hor high explosives). Let us 
now specialize these results to equatorial Spread F bubbles. In this 
case, neglecting for the moment any horizontal ambient electric field, 

A 

E in Fig. 1 is really (ftM/ev^) £x b and it is the polarization 
(induced) electric field (in the direction of E_' and equal to (flM/ev. ) 

A 

(6n/(n - <5n)) £x b which causes the bubble to rise with respect to 

the ionosphere via the local x convective drift motion. The picture 
exhibited in Fig. 1 for Spread F bubbles (using essentially eq. (4) 
which contains collisional Rayleigh-Taylor effects) shows that the 
bubbles spread more vertically than horizontally. This agrees with the 
radar observations [ Woodman and La Hoz , 1976] and of course the numerical 
simulations [ Scannapieco and Ossakow , 1976] which utilized eq. (1) and (2). 
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Linear Spread F Bubble Rise Velocity Model 




Let us now look at Spread F bubble rise velocities for various 

bubble models. We will utilize eq. (4) and in our first case let us 

include a horizontal ambient electric field, E , to show its effects. 

—o 

Then the equation becomes 


V (nV <J>) = E* • V n 

i' i ' i 



E*eE+ — - jxb 
-° v in e 


Our first model will be a linearized model such that 


n = n (x) + ni e'(ky-wt) 





i(ky-wt) 


where n (x) is the equilibrium altitude dependent density profile and 
subscript one refers to perturbed quantities. Substituting (8) into (7) 
and taking E to be in the y direction yields for the polarization 
potential 


. = * . fl M , r \ i n i 

^ 1 " 1 (v • e ^ o) kn 
in ' o 


The bubble velocity is 
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(10) 



In eq. (10) we use £ = —V4>^ = - i]^ so that for the linear case |Vg 
|cE/B| = |ck<|)j/Bj and we have for the vertical velocity (6n e n^) 



( 11 ) 


(Note: Here 6n is positive and taken to mean the magnitude of the 

depletion.) It is interesting to note that adding E in the direction 

0 


parallel (antiparallel) to £ x b increases (decreases) the polarization 

field (as given by ) and makes the bubbles rise faster (slower). The 

numerical simulations of Scannapieco and Ossakow L1976] did not have a 

horizontal electric field and adding one in the proper direction would 

make the bubble rise faster. However let us look at the effect of E . 

We wish to compare the two terms inside the parenthesis of eq. (11). If 

we assume E 'v* lmV/m, B = 0.3 gauss then cE Q /B 'v 33 m/sec. At 350 km 

altitude g/v. ^ 40 m/sec, whereas at 450 km g/v.- n 'v 220 m/sec. Thus 

we see that as one goes to higher altitudes E makes a more negligible 

contribution to the bubble rise velocity. Indeed, even if the electric 

field reverses, once the bubble gets to higher altitudes (decreasing 

v. ) the bubble would keep rising. This bubble rising even after the 
in 

electric field reverses has been seen experimentally [Kelley, private 
communication, 1976]. 
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(12) 


Neglecting E in eq. (11) we have 


Using Fig. 2 for v. as a function of altitude, Fig. 3 depicts V« as a 
function of altitude for various values of percentage depletions, 6n/n Q 
(large 6n/n represents usage of eq. (12) beyond the range of validity 
and perhaps represents a lower bound). Immediately it can be seen from 
this figure that large bubble velocities (^100 m/sec) requires large 
depletions and/or high altitudes. 

Nonlinear Spread F Bubble Rise Velocity Model 

In the case of a so-called sheet geometry in which we allow n and 
4> to be functions of y only, eq. (4) can be integrated once to yield 


ay 



q K 
v in e 


(13) 


where n is the ambient density outside the bubble. For this case the 
0 

dimension of the bubble in the vertical direction is extremely long 
compared with its horizontal dimension. Using eq. (13), the vertical 
bubble rise velocity now becomes 




, 6n 
n 


v. 
l n 


o 


(14) 
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where we have used Vg = 1^-y x b and n = n Q - 6n for a depletion. A 

peculiarity of this model is that V g °° as 6n/n 1, i.e.; for a 100% 
depletion. Fig. 4 shows the bubble rise velocity as a function of alti¬ 
tude for various vales of 6n/n , using eq. (14) and values of v. depicted 

o in 

in Fig. 2. Here, as in Fig. 3, high altitudes and/or large depletions 
yield high vertical bubble velocities. 

As we have stated previously, eq. (3) and (4) in general do not 
admit of a two dimensional (i B) solution. However, there is a patho¬ 
logical case where one can obtain exact solutions and we will discuss 
this now. We write for the density 


n(x,y) = n Q 


n D (x,y) = n Q 



(15) 


r 


where n is the ambient ionospheric density, n g represents the depletion 
or bubble (6n/n is the magnitude of the ratio of maximum depleted density 
to ambient density), and F describes the shape of depletion. If n g (x,y) 

•+ 0 as x,y °° and V n is finite then there is no equilibrium solution 
to eq. (3) and (4). However, if the depletion density is a constant and 
the shape of the density depletion contour is an ellipse described by 



[where H(x) is the Heaviside step function and a and b represent the 
vertical dimension (axis of the ellipse) and horizontal (east-west) 


10 









dimension respectively], then the electric field inside the depletion 
is a constant and the depletion (bubble) moves with constant velocity 
without changing its shape. This method was first employed by Linson 


[1972] for barium clouds (plasma density enhancements) where eq. (5) and 
(6) were utilized. (Note: In this case the minus signs in eq. (15) 
would become plus signs.) Indeed, eq. (4) or (6) can be thought of as 
the case of a dielectric immersed in a uniform electric field. However, 
the plasma depletion is analogous to the case of a cavity immersed in a 
dielectric with a uniform field (since the field is enhanced inside a 
depletion). The solution of the potential eq. (4) or (6) for the general 
elliptical shape can be found in Smythe [1950]. For the case of a plasma 


depletion described by eq. (4), with £ x b in the y direction, we have 
for the polarization (induced) electric field 



(16) 


with the attendant bubble rise velocity given by 



(17) 
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(Note: In arriving at eq. (16) we have assumed v. to be a constant. 
This is a good approximation so long as the size of the depletion is 
small compared with an atmospheric scale height). 

Now we will examine eq. (17) with respect to different shapes. For 
the case b « a, i.e., the sheet geometry, we obtain 


V 


B 





which is eq. (14) as it should be. For a/b ->-0, i.e., a slab case 

A 

(uniform in the direction of £ x b), V g -+■ 0 since the induced electric 
field goes to zero. When a = b (cylindrical case), we have 


V 


B 


6n 




(18) 


For 5n/n 0 « 1 eq. (18) yields V g = 1/2 (g/v^ )(6n/n ) which is half 

the value given by the linear perturbation model (see eq. (12)). With 

Sn/n = 1 (100% depletion), eq. (18) yields V g = g/v^. Fig. 5 shows 

the vertical bubble rise velocity for a cylindrical shape (as given by 

eq. (18)) and Fig. 6 exhibits the bubble rise velocity for an elliptical 

shape (a/b = 5) as a function of altitude for various values of 6n/n , 

using values of v. depicted in Fig. 2. 

3 in 
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III. DISCUSSION AND SUMMARY 


Figures 3-6 depict vertical bubble (plasma depletion) rise veloci¬ 
ties for various bubble models (shapes), for equatorial Spread F geo¬ 
metry, as a function of altitude and percentage depletions. Exceedingly 
high velocities can be attained, especially for the severely elongated 
models (e.g., sheet model in Fig. 4 and the elliptical model with a 5:1 
axial ratio in Fig. 6). Table 1 shows the bubble rise velocity, in 
units of g/v in , as a function of fractional depleted plasma density for 
various bubble shapes. For all values of 6n/n Q (except at 1.0) the non¬ 
linear cylindrical model gives a smaller velocity (see also eqs. (12) 
and (18)) than the linear model. The table also shows that the cylin¬ 
drical bubble velocity for a given 6n/n Q is attained by the sheet model 
for a value of 1/2 6n/n (see also eqs. (14) and (18)). An elliptical 
model with 10:1 axial ratio is also shown in Table 1. Except for 
6n/n = 1, it has already become fairly sheet-like. A common feature to 
all models is that the higher the altitude and/or the higher the density 
depletion, the higher the vertical rise velocity. This result is con¬ 
sistent with the recent satellite in situ measurements described by 
McClure et al. [1977]. 

The models described by Fig. 3-6 essentially assume that the bubble 
does not change its shape. For real bubbles, where the depletion den¬ 
sity changes inside the bubble, the bubble should elongate in the 
vertical direction (steepening on the topside of the bubble) and become 
narrower in the east-west extent. This is in agreement with the radar 
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observations ( Woodman and La Hoz , 1976) and the numerical simulation of 
rising equatorial Spread F bubbles ( Scannapieco and Ossakow , 1976). 
Therefore, in modeling a changing bubble by one of constant shape we are 
employing somewhat of an artifact. However, using this technique Linson 
[1972] has shown (using optical data) that evolving barium clouds (plasma 
density enhancements) tend to be more cylindrical than sheet-like. Using 
the constant shape approach and the simulation results of Scannapieco 
and Ossakow [1976], where V~ ^ 10 m/sec and 6n/n Q ^ 40% for an altitude 
of 350 km, the cylindrical bubble shape is favored (although the simula¬ 
tion pictures show an elongated bubble). More data and simulation 
results are needed to ascertain which model fits best. 

All of the models described in this paper, based on collisional 
Rayleigh-Taylor type equations, show that the vertical bubble rise 
velocity is of the form 



(19) 


where f(6n/n ) is some function of 6n/n Q . The factor g/v.j n should not 
be too surprising since the growth rate (y) for the collisional Rayleigh- 
Taylor instability, neglecting recombination effects, is g/v. L ( Hudson 
and Kennel , 1975; Scannapieco and Ossakow , 1976), which is a velocity 
divided by the electron density gradient scale length. This of course 
is similar to the growth rate (y = cE /BL) for barium cloud striations 












due to E„ x B gradient drift instability (Simon, 1963; Linson and Workman, 

O - " - - ■ ■ • 

1970). This similarity between the gravitational and the electric field 
terms has been shown in eq. (11). In eq. (19), the factor g/v^ favors 
high altitude because v. decreases as a function of altitude (see 
Fig. 2). The factor f increases with increasing 6n/n Q (see Fig. 3-6). 

At least according to linear theory, 6n/n Q increases with altitude and 
decreasing bottomside electron density gradient scale length (y = g/v^ L). 
These ideas indicate that there is a range of bubble vertical rise 
velocities (in accord with McClure et al. , [1977]), but these rise 
velocities depend on ambient equatorial F region ionospheric conditions 
(e.g., height of the F peak and bottomside gradient scale length). 

Recent numerical simulations (S. Ossakow, S. Zalesak and B. E. McDonald, 
private communication, 1977) support these ideas. Also, it should be 
pointed out that for 6n/n = 0, i.e., for bubbles which have the same 
density as the background ionosphere, the induced (polarization) electric 
field (which causes the bubble to rise through E x B) becomes zero (with 
all models, as it should), as does the attendant vertical bubble rise 
velocity. This simply states that a bubble created on the bottomside of 
the F region will rise until its density equals the topside background 
ionospheric density. 

It should be pointed out that our vertical velocity result for 
cylindrical bubbles does not agree with Ott's [1977] result. In his 
analogy with bubbles in a fluid, Ott found only a cylindrical solution. 

For the collision dominated regime he found Vg = l/2(g/v. )(6n/n ), in 
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our terminology; whereas, our result shows Vg = 1/2(g/v^. n )(6n/n Q ) • 

(1 - 1/2(6n/n Q )) . For the case where the bubble density is negligibly 

small (6n/n ~ 1) compared with the density in the ionosphere outside, 
Ott's result yields Vg = l/2(g/v. ); whereas, our result yields 

V D = g/v. . 

B in 

Finally some comments about the various upward and westward bubble 
velocities presented by McClure et al. [1977] and the tilts of the 
plumes seen by the Jicamarca radar measurements ( Woodman and La Hoz , 

1976). Equations (1) and (2), used by Scannapieco and Ossakow [1976], 
or eqs. (3) and (4) for that matter produce only vertical bubble rise. 
These equations only admit of a horizontal polarization electric field 

A 

(eastward) caused by the cj. x b term such that depletions rise and enhance¬ 
ments fall with respect to the ionosphere. Inclusion of an ambient 
background horizontal electric field (eastward or westward) in the equa¬ 
tions could only cause additional vertical (upward or downward) motion 
of the bubbles. Since the ionosphere moves mainly downward at night 
(westward electric field; see McClure et al. , 1977) this could slow down 
the bubbles at low altitudes and have a negligible effect at higher 
altitudes (see eq. (11)). Vertically the bubbles move more slowly at 
lower altitudes and/or small values of 6n/n . Thus there exist conditions 
whereby the bubbles do not move vertically with respect to the iono¬ 
sphere. McClure et al. [1977] and Woodman and La Hoz [1976] have 
stated that the ionosphere moves mainly eastward at night. This would 
mean an ambient electric field pointing downward. Inclusion of such an 
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electric field in the equations would cause bubbles to move eastward. 

Thus the westward motion of the bubbles must be caused by introducing 
an asymmetry into the eq. (1) and (2) (or (3) and (4)) which causes a 
vertical (upward) electric field. Such an effect can be produced by an 
eastward directed neutral wind which would essentially polarize the 
bubble to produce an upward polarization electric field and account for 
the westward motion of the bubble (and tilts of the radar plumes). This 
last idea has also been set forth by Woodman and La Hoz [1976] and Ott 
[1977]. 

Recalling our discussion about the electric field, E* = E„ + (ft/v. ). 

-o' in' 

A 

(M/e) £ x b, after eq. (6), all of our vertical bubble rise velocities can 
take account of horizontal motions. The general bubble velocity is then 
just Vg = (c/B) E* x b f(<Sn/n Q ) = ((c/B) x b - £/v in ) f(6n/n Q ). This 
clearly shows how the directionality of Vn is formed. In this case for 

A 

the elliptical shape a is in the direction of E* x b, etc. and non¬ 
constant density depletions will steepen and elongate in this direction. 

Recent results from the Wideband associated equatorial experiments 
performed at Kwajalein (August 1977) lend credence to this notion ( $. 
Ossakow and J. Clark , private communication, 1977). Backscatter irre¬ 
gularity ('vl meter) measurements (performed by D. Towle) have shown that 
these patches rise and appear westward when neutral wind measurements 
(Fabry-Perot interferometer measurements performed by M. Biondi and D. Sipler) show 

strong west to east winds (M60 m/sec) earlier in the evening. As the 
evening progresses the neutral winds die down (M0-20 m/sec) and the 
radar shows the irregularities appearing eastward. 
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Table 1. Bubble rise velocity (in units of g/v. ), VA, as a function 
- in d 

of fractional depleted density, 6n/n Q , for various bubble shapes. 
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(A) 

PLASMA DENSITY 
ENHANCEMENT 


(B) 

PLASMA DENSITY 
DEPLETION 


(a) t «0 





Fig. 1 — Deformation of (A) a plasma density enhancement (e.g., a barium 
cloud) and (B) a plasma density depletion (e.g., equatorial Spread F bubble 
or artificially created “hole”) in crossed electric (E^) and magnetic (B) fields 
(representative of equatorial geometry), (a) Initial (t=0) cylindrical isoden¬ 
sity contours with largest enhancement (+) and depletion (-) values in the 
center and falling off away from these values to an ambient level n Q outside. 
E/ represents the total electric field inside the enhancement or depletion 
and is the ambient field outside, (b) Qualitative evolution of the isoden¬ 
sity contours, (c) Approximate values for the initial electric field (E') inside 
the enhancement or depletion. 


23 








ALTITUDE (km) 



Fig. 2 — Ion-neutral collision frequency, as a 
function of altitude 
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Fig. 3 — Depletion (bubble) vertical rise velocity, V B , 
as a function of altitude for various values of the per¬ 
centage depletion, 5n/n 0 , for the linear bubble model, 
i.e., V B = (g/u in )(6n/n 0 ) 
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Fig. 4 — Same as fig. 3, but for the nonlinear sheet bubble 
model, i.e., V B = (g/u in )(6n/n 0 )(l-6n/n 0 )- 1 
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Fig. 5 — Same as fig. 3, but for the nonlinear cylindrical 
bubble model, i.e., V B = (g/u in )(5n/n 0 )(2-5n/n 0 )' 1 
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Fig. 6 — Same as fig. 3, but for the nonlinear elliptical 
bubble model (with the vertical dimension five times the 
horizontal dimension), i.e., V B = (g/u in )(5 5n/n 0 )(6-5 
5n/n 0 )'i 
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